Double barrier potentials for matter-wave gap solitons 



OO 
O 

o 






a 



I 

o 
o 



> 

00 
(N 

m 

^* 

o 

00 
O 

> 

X 

S-H 



V. Ahufinger/ B. A. Malomed,^ G. Birkl,^ R. Corbalan,^ and A. Sanpera^ 

^ ICREA and Grup d' Optica, Departament de Fisica, 

Universitat Autdnoma de Barcelona, E-08193 Bellaterra, Barcelona, Spam. 

'^Department of Physical Electronics, School of Electrical Engineering, 

Faculty of Engineering, Tel Aviv University, Tel Aviv 69979, Israel 

^ Institut fiir Angewandte Physik, Technische Universitat Darmstadt, 64289 Darmstadt, Germany 

''Grup d'Optica, Departament de Fisica, Universitat Autdnoma de Barcelona, E-08193 Bellaterra, Barcelona, Spain. 

^ICREA and Grup de Fisica Tedrica, Departament de Fisica, 

Universitat Autdnoma de Barcelona, E-08193 Bellaterra, Barcelona, Spam. 

We investigate collisions of solitons of the gap type, supported by a lattice potential in repulsive 
Bose-Einstein condensates, with an effective double-barrier potential that resembles a Fabry-Perot 
cavity. We identify conditions under which the trapping of the entire incident soliton in the cavity 
is possible. Collisions of the incident soliton with an earlier trapped one are considered too. In the 
latter case, many outcomes of the collisions are identified, including merging, release of the trapped 
soliton with or without being replaced by the incoming one, and trapping of both solitons. 



I. INTRODUCTION 

Bright matter-wave solitons [2 |2] provide an excep- 
tional testbed for studying quantum mechanics above the 
single-atom level. For the observation of quantum phe- 
nomena, coherence is a crucial requirement, and matter- 
wave solitons, which propagate without dispersion, may 
render this observation more feasible. 

The generation of solitary waves in quasi-one- 
dimensional (QID) attractive Bose-Einstein condensates 
(BEG) [1 may be regarded as a tuned equilibrium be- 
tween the dispersive effects that tend to spread the 
atomic wave function, and the nonlinear attractive inter- 
actions which oppose the spreading by providing an effec- 
tive self-focusing of the matter waves. The consequence 
of such an equilibrium is a stable mesoscopic atomic wave 
packet propagating without dispersion. 

Repulsive condensates in QID geometries can support 
matter- wave solitons of the gap type j2] if they are loaded 
in an optical lattice (OL), which gives rise to the bandgap 
structure, and are placed at the edge of the first Brillouin 
zone, where there is a gap between the first and second 
bands \W . Under these conditions, the soliton exhibits 
an effective negative mass, which permits to balance the 
dispersion and the nonlinear interactions, even if they 
are repulsive. Various effects generated by the negative 
effective mass were considered in Refs. [21 [5l El |7] . 

Stabihty conditions for gap solitons (GSs) impose se- 
vere restrictions on their interactions with a potential 
well or barrier corresponding to a local modification of 
the periodic structure. Specifically, the requirement of 
the stability allows only for perfect transmission or per- 
fect reflection, but not partial transmission and reflection 
that plane waves display [5]. In other words, GSs can- 
not split through the interaction with linear defects, and 
behave like particles exhibiting mesoscopic quantum fea- 
tures [6 , including the quantum reflection of the entire 
soliton - an effect that has been also reported recently for 
matter- wave solitons in the self-attractive BEG case [9]. 
Although quantum reflection of ultracold atoms from a 



solid surface has been reported too [TU] , the limit of the 
complete (100%) reflection, predicted for the solitons, is 
not achievable in that case. Interactions of matter-wave 
solitons with nonlinear traps and barriers produced by 
spatial variations of the scattering length have also been 
recently addressed [TT| . 

In this work we aim to explore such quantum features 
in the case of the interaction of a matter-wave soliton of 
the gap type with a double-barrier potential resembling 
a Fabry-Perot cavity. In particular, we demonstrate that 
it is possible to trap a soliton in the cavity formed by the 
two potential barriers (a similar effect for optical solitons 
in fiber Bragg gratings was predicted in Ref. ilJj). In 
the context of BEG, the propagation through a double- 
barrier potential acting as a Fabry-Perot interferometer 
for matter-waves leads to bistability of the transmitted 
flux and resonant transport |13j . 

The paper is organized as follows. In Sec. [IT] we formu- 
late the physical model which includes the effective cavity 
for the GS. Sec. |III| is devoted to the study of the condi- 
tions for the trapping of the entire soliton in the cavity. 
In Sec. |IV| the collision of a second GS with one trapped 
in the cavity is studied (this interaction also bears some 
similarity of collisions between free and defect-trapped 
optical solitons in fiber gratings [Mj). The paper is con- 
cluded in Sec. El 



II. THE PHYSICAL MODEL 

The dynamics of bright GSs created in a QID geometry 
at zero temperature may be accurately described by the 
one-dimensional (ID) Gross-Pitaevskii equation (GPE): 



dt 



- — A + F(x)H-ff|V|' 



^, 



(1) 



where the effective nonlinearity is g = 2hasLiJt , with Og the 
s-wave scattering length and ujt the transverse trapping 



frequency. The effective axial potential, 

V [x] = (l/2)mw2^2 _^ y^ sin^{-Kx/d), (2) 

includes the parabolic trap with corresponding frequency 
ojj, and the OL with spatial period d and depth Vq- 

First, we briefly summarize the numerical procedure 
used to generate GSs, as per Ref. [5]. The ground 
state is found for a *^Rb condensate (a^ = 5.8 nm) 
formed by A^ = 500 atoms trapped magnetically with 
transverse and axial frequencies oJt = 715 x 27r Hz and 
Wa; = 14 X 27r Hz, respectively, in the presence of an 
OL, with potential depth equal to the OL recoil energy, 
Vq = Er = h^k^ /2m, where k = n/d is the recoil mo- 
mentum, and d = 397. 5nm, the period. Then, the axial 
magnetic trap is suddenly turned off and an appropri- 
ate phase imprinting leads to the inversion of the sign of 
the wave function at each second site. As a result, the 
system evolves toward a self-sustained staggered soliton 
with a negative effective mass. The soliton, which con- 
tains approximately 35% of the initial number of atoms 
and extends over ~ 11 sites of the OL potential, is gen- 
erated at rest with the center at x = 0. In order to set it 
into motion, we instantaneously impart an appropriate 
momentum to the soliton, to which it responds by self- 
adaptation to the new conditions, i.e., expunging atoms 
until a new equilibrium state is reached. With momen- 
tum p = 0.1kh lent to the GS, it settles down into a state 
with 27% of the initial number of atoms (-/Vgnai — 135), 
total energy 0.92£'r, and the kinetic energy of its motion 
Ek = QMEr [6]. 

Our aim is to study the interaction of the so generated 
moving matter-wave soliton with a Fabry-Perot type po- 
tential, formed by two potential barriers forming a cavity 
for the soliton, cf. a similar configuration proposed for 
optical solitons in fiber gratings 112] . After turning off 
the magnetic trap, V{x) in Eq. M accounts only for the 
OL, which is locally modified around two sites, Xmi and 
2^7712, to generate the double-barrier structure as follows: 



V {x) = Vq sin^^TTx/d) + V;nod (x) . 



(3) 



where Vmod (a;) reads: 



(x- 




{x- 



if a^mi 

if a;m2 
otherwise 



l/2<x< x^i + 1/2; 

1/2<X< Xm2 + 1/2] 



(4) 

with a — 6d. Below, we consider the symmetric cavity 
created by two identical barriers, with Vmi — Vm2 — 
Vm < 0. Note that, due to the negative effective mass 
of the GS, the barriers actually correspond to a local 
decrease of the periodic potential. Points Xmi,2 are fixed 
at local minima of the OL potential, and the distance 
between the barrier centers, A = Xm2 — Xmi, is assumed 
large enough to have the actual size of the cavity, B = 
A — I, much larger than the axial size of the soliton. 
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FIG. 1: Diagram of outcomes of the collisions of a moving 
lattice soliton with a cavity formed by two identical barri- 
ers: reflection, transmission, and trapping, in the white, light 
grey, and dark grey areas, respectively. The effective size of 
the cavity is B = 20d. Parameters / (measured in units of 
the lattice period, d) and Vm, measured in units of the recoil 
energy, E,-, are the width and height of the two barriers, re- 
spectively. The horizontal dotted line designates the kinetic 
energy of the moving soliton. 



III. THE CAVITY 

The interaction of a moving matter-wave GS with a 
single barrier in the OL was addressed recently in Ref. 
f6], where it was shown that the soliton does not split. 
For a fixed kinetic energy, there exists an abrupt transi- 
tion from complete transmission to complete reflection, 
as the height of the barrier increases, for all considered 
values of the barrier widths. The border between these 
two outcomes of the collision of the soliton with the 
barrier is shown, in the plane of the barrier's parame- 
ters, height Vm and width I, by filled circles in Fig. [T] 
Complete reflection and transmission occur, respectively, 
above and below this border (i.e., the white area in the 
figure corresponds to the bounce of the entire soliton from 
the barrier). In the process of the collision, the soliton 
naturally decreases its velocity (even if does not bounce). 
This slowing down of the soliton in the region of the de- 
fect is more pronounced as we approach the border be- 
tween the two behaviors. Addressing the configuration 
with the second barrier placed at a certain distance from 
the first one, we notice that the behavior of the incident 
soliton which hits the first barrier remains as described 
before, i.e., a sudden transition from perfect transmis- 
sion to perfect refiection occurs. Nevertheless, when the 
height of the barriers is close to the transition point, and 
the length of the cavity (B) is larger than the size of the 
soliton, the soliton which has passed the first barrier gets 
trapped in the cavity, in a state of oscillatory motion. 
Thus, three scenarios are identified in the interaction 




FIG. 2: Contour plots of the spatiotemporal evolution of the 
lattice soliton colliding with the set of two identical barriers 
centered at Xmi = 25d and Xm2 ~ 45d, each with width 
I = 2d and the following height: (a) \Vm\ = O.OOSBr; (b) 
\Vm\ = O.Oll^r; (c) \Vm\ = O.OUEr. Here as well as in the 
following figures the horizontal dotted lines show the position 
of the centers of both barriers. 



of the incident soliton with the double-barrier structure: 
complete reflection, complete transmission, and trapping 
into the oscillatory state. As said above, the reflection 
occurs in the same area of the parameter space (white 
region in Fig. [l]) as for the single barrier. The com- 
plete transmission takes place for values of the barrier's 
height and width well inside the transmission region for 
the single barrier (the light grey area in Fig. IT]) , while the 
trapping is observed close to the transmission-reflection 
border for the single barrier (the dark grey region in Fig. 
fl]). The results shown in Fig. fl] correspond to a flxed 
cavity size of B = 20d and we have checked that the 
parameter space for trapping slightly increases with the 
size of the cavity. For a fixed width of the barriers, and a 
large enough distance between them, the three scenarios 
follow each other with the increase of the barriers' height. 
These scenarios are illustrated, in Fig. |2]by spatiotem- 
poral trajectories of the soliton hitting the cavity formed 
by two identical barriers of width I = 2d, which are sepa- 
rated by distance A — 20d, giving the cavity enough room 
to trap the soliton, B — I8d. The complete transmission 
is displayed in panel (a) for \Vm\ = O.OOSEr, panel (b) 
with \Vm\ = O.Ollii^r shows an example of the trapping, 
and the bounce (complete reflection) is observed in panel 
(c), for \Vm\ = 0.014i^^r. Figure|2jb) clearly demonstrate 
that the trapped soliton performs periodic oscillations in 
the cavity without any visible loss of atoms. The period 
of the oscillations, for given parameters of the barriers, 
can be modified by changing the size of the cavity, as 
shown in Fig. [3) Note that the smallest period, lim- 
ited by the condition that the size of the cavity must be 
larger than the soliton's axial size, is ~ 30 ms, in physi- 
cal units. Faster oscillations were predicted for a lattice 
soliton trapped in a potential well [6 . It is relevant to 
note that, in the above-mentioned case of the trapping 
of optical solitons by a cavity in the model of the fiber 
grating |12j . the maximum capture efficiency (in terms 
of the soliton's energy) is no more than 60%, contrary to 
what happens with the matter-wave GSs, which may be 
trapped entirely, without losses. 
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FIG. 3: Examples of the trapping of the lattice soliton by the 
pair of barriers of width I = 2d and height \Vm\ = 0.011-Er. 
The first barrier is centered at Xmi = 25d, and the second at; 
(a) Xm2 — 45d, {h)xm2 = 55d, and (c) Xm2 — 90d. 



IV. COLLISIONS BETWEEN FREE AND 
TRAPPED SOLITONS 



The next natural step in the analysis is to consider a 
collision between an incident soliton with the cavity al- 
ready occupied by an (identical) earlier trapped GS. Dif- 
ferent outcomes of the collisions are observed, depending 
on time delay At between the two solitons. Figure |4] dis- 
plays three cases for the cavity formed by two barriers of 
width I — 2d and height \Vm\ = O.OllSr, separated by a 
distance oi A — 20d, corresponding to (a) At — 20 ms, 
(b) At = 30 ms, and (c) At — 35 ms. In (a), the incident 
soliton bounces back, while the trapped one performs os- 
cillations in the cavity; in (b), the two solitons merge 
into a single one, and in (c), the incident soliton bounces 
back, kicking out the trapped one in the forward direc- 
tion. Two more cases are shown in Fig. |4]for the same 
parameters of the barriers, but for a different separation 
between them, A = 30d, and time delays i = 25 ms in 
(d), and i = 30 ms (e). In this case, the cavity has enough 
room to trap the two solitons, which gives rise to new col- 
lision scenarios: in (d), the incoming soliton gets trapped 
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FIG. 4: Different outcomes of the collision between an inci- 
dent GS with an identical one previously trapped. The barri- 
ers of width I = 2d and height \Vm\ = 0.011i?r are separated 
by distance A = 20d, in (a), (b) and (c), or yl = 30d, in (d) 
and (e). The second soliton is send after the first one with 
time delay At = 20 ms (a), At = 30 ms (b), At = 35 ms (c), 
At = 25 ms (d), and At = 30 ms (e). 



by kicking out the previously trapped one ("recharge"), 
while in (e) both solitons get trapped in the cavity, oscil- 
lating in counter-phase. These sundry dynamical behav- 
iors suggest new experimental possibilities for the control 
and manipulation of matter-wave GSs. 

V. CONCLUSIONS 



In this work, we have studied the interactions of bright 
matter-wave solitons of the gap type, with negative ef- 
fective mass, which are supported by the interplay of the 
OL (optical lattice) and repulsive nonlinearity in BEG, 
with a cavity formed by two far separated identical local 
potential barriers. We have shown that there exists a 
parameter region in which the incident soliton is trapped 
by the cavity into the shuttle state. This region can be 
found for all the values of the barriers' width, provided 
that their height corresponds to the transmission of the 
soliton by a single barrier, but close to the reflection- 
transmission border. 

The interaction of a second soliton which hits the cav- 
ity already occupied by a trapped oscillating soliton has 
been considered too. In that case, a number of differ- 
ent collision scenarios can be identified, depending on 
the time delay between the launch of the two solitons. 
Particularly interesting outcomes are the merging of the 
two solitons at the position of the first barrier, bounce 
of the second soliton kicking out the trapped one, the 
"recharge", i.e., release of the originally trapped soliton 
which is replaced by the incident one, and trapping of 
both solitons into the state of shuttle oscillations in the 
cavity, with a phase shift of tt between them. 
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